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Abstract. We have calculated diffractive dijet production in deep-inelastic scattering (DIS) at low- 
Q 2 and next-to-leading order (NLO) of perturbative QCD, including contributions from direct and 
resolved photons. We study how the cross section depends on the factorization scheme and scale My 
at the virtual photon vertex for the occurance of factorization breaking. The strong My-dependence, 
which is present when only the resolved cross section is suppressed, is tamed by intodrucing the 
suppression also in the initial-state NLO correction of the direct part. 



INTRODUCTION 

From a perturbative QCD (pQCD) point of view, the central question for hard diffractive 
scattering events, characterized by a large rapidity gap devoid of particles in high-energy 
collisions, is whether they can be factorized into non-perturbative diffractive parton 
density functions (PDFs) of a colorless object (e.g. a pomeron) and perturbative partonic 
cross sections. This concept, based on a long-standing proposal by Ingelman and Schlein 
J3], is believed to hold for the scattering of point-like electromagnetic probes off a 
hadronic target, such as deep-inelastic scattering (DIS) or direct photoproduction J2], 
whereas it has been shown to fail for purely hadronic collisions 

Factorization is thus expected to fail also in resolved photoproduction, where the 
photon first dissolves into partonic constituents, before these scatter off the hadronic 
target. The separation of these two types of photoproduction events is, however, a 
leading order (LO) concept. At next-to-leading order (NLO) of pQCD, they are closely 
connected by an initial state singularity originating from the splitting y — > qq 0], 
which may play a crucial role in the way factorization breaks down in diffractive 
photoproduction [5]. Factorization breaking effects are therefore expected to show up 
first in observables that distinguish between direct and resolved photoproduction, such 
as distributions in the longitudinal momentum fraction Xy of partons in the photon JfJ, 
the photon virtuality Q 2 CD, or the dependence of the predicted cross sections on the 
factorization scale My J5J. It is clear that this M r -dependence is unphysical and must 
be remedied also for the case of factorization breaking of the resolved part of the cross 
section. A proposal how to achieve this has been worked out in our previous work ||5j] and 
will be described in the next Section. For demonstrative purposes we restrict ourselves 
to low-<2 2 DIS using the kinematic framework of our earlier publication CD. 



Three groups have recently tried to extract diffractive parton densities from inclusive 
diffractive DIS data at DESY HERA, treating the Pomeron either as a hadronic object 
within Regge factorization JalSI] 

fa(x,Q 2 ;XP,t)=fp/ p (xp,t)f a /p(P=x/xp,Q 2 ) (1) 

or perturbative QCD [10], but so far only the HI set has been tested for factorization in 
a different scattering environment, i.e. in photoproduction of dijets at HERA 111 ill . 



FACTORIZATION SCHEME AND SCALE DEPENDENCE 

A factorization scheme for virtual photoproduction has been defined and the full NLO 
corrections for inclusive dijet production have been calculated in [ 12]. They have been 
implemented in the NLO Monte Carlo program JET VIP II 1 311 . We have adapted this NLO 
framework to diffractive dijet production. According to Il2ll . the subtraction term, which 
is absorbed into the PDFs of the virtual photon f a /y(xy,My), is of the form as given in 
|@]. The main term is proportional to \n(M 2 /Q 2 ) times the splitting function 

7 2_i_fi __\2 

P qi ^y(z)=2N c Q f +l 2 Z) , (2) 

where z = PiPi/poq £ [*', 1] and Qi is the fractional charge of the quark qi. p\ and p2 
are the momenta of the two outgoing jets, and po and q are the momenta of the ingoing 
parton and virtual photon, respectively. Since Q 2 = —q 2 <C Mi, the subtraction term 
is large and is therefore resummed by the DGLAP evolution equations for the virtual 
photon PDFs. After this subtraction, the finite term M(Q 2 )y^, which remains in the 
matrix element for the NLO correction to the direct process Ill2ll . has the same My- 
dependence as the subtraction term, i.e. InMy is multiplied with the same factor. As 
already mentioned, this yields the My-dependence before the evolution is turned on. In 
the usual non-diffractive dijet photoproduction these two My-dependences cancel, when 



the NLO correction to the direct part is added to the LO resolved cross section II 1411 . 
Then it is obvious that the approximate M r -independence is destroyed, if the resolved 
cross section is multiplied by a suppression factor R to account for the factorization 
breaking in the experimental data. To remedy this deficiency, we propose to multiply the 
In My-dependent term in M(Q 2 )y^ with the same suppression factor as the resolved cross 
section. This is done in the following way: We split M(Q 2 )y^ into two terms using the 
scale pj in such a way that the term containing the slicing parameter y s , which was used 
to separate the initial-state singular contribution, remains unsuppressed. In particular, 
we replace the finite term after the subtraction by 



M(Q 2 ,R] 



MS 




(3) 



where R is the suppression factor. This expression coincides with the finite term after 
subtraction (see Ref. [5]) for R = 1, as it should, and leaves the second term in Eq. 
© unsuppressed. In Eq. © we have suppressed in addition to \n(M 2 /p^) also the 

z-dependent term ln(z/(l — z)), which is specific to the MS subtraction scheme as 
defined in 1 1211 - The second term in Eq. © must be left in its original form, i.e. 
being unsuppressed, in order to achieve the cancellation of the slicing parameter (y s ) 
dependence of the complete NLO correction in the limit of very small Q 2 or equivalently 
very large s. It is clear that the suppression of this part of the NLO correction to the direct 
cross section will change the full cross section only very little as long as we choose 
My ~ pj. The first term in Eq. ©, which has the suppression factor R, will be denoted 
by DIRis in the following. 

To study the left-over My-dependence of the physical cross section, we have calculated 
the diffractive dijet cross section with the same kinematic constraints as in the HI 
experiment 11511 . Jets are defined by the CDF cone algorithm with jet radius equal to 
one and asymmetric cuts for the transverse momenta of the two jets required for infrared 
stable comparisons with the NLO calculations 111 611 . The original HI analysis actually 
used a symmetric cut of 4 GeV on the transverse momenta of both jets [ 17]. The data 
have, however, been reanalyzed for asymmetric cuts IU51 . 

For the NLO resolved virtual photon predictions, we have used the PDFs SaSID 
iflill and transformed them from the DIS r to the MS scheme as in Ref. lfl2ll . If not 
stated otherwise, the renormalization and factorization scales at the pomeron and the 
photon vertex are equal and fixed to p\ = p* T - etV We include four flavors, i.e. rif = 4 
in the formula for a s and in the PDFs of the pomeron and the photon. With these 
assumptions we have calculated the same cross section as in our previous work Q]. First 
we investigated how the cross section da/dQ 2 depends on the factorization scheme of 
the PDFs for the virtual photon, i.e. do/dQ 2 is calculated for the choice SaSID and 
SaSlM. Here da/dQ 2 is the full cross section (sum of direct and resolved) integrated 
over the momentum and rapidity ranges as in the HI analysis. The results, shown in 
Fig. 2 of Ref. [5] demonstrate that the choice of the factorization scheme of the virtual 
photon PDFs has negligible influence on da/dQ 2 for all considered Q 2 . The predictions 
agree reasonably well with the preliminary HI data 11511 . 

We now turn to the My-dependence of the cross section with a suppression factor 
for DIRis, which is the main part of this Report. To show this dependence for the two 
suppression mechanisms, (i) suppression of the resolved cross section only and (ii) ad- 
ditional suppression of the DIRis term as defined in Eq. © in the NLO correction of 
the direct cross section, we consider da/dQ 2 for the lowest <2 2 -bin, Q 2 6 [4,6] GeV 2 . 
In the left part of Fig. [H this cross section is plotted as a function of £ = My/ pj in the 
range £ £ [0.25; 4] for the cases (i) (light full curve) and (ii) (full curve). We see that 
the cross section for case (i) has an appreciable t, -dependence in the considered t, range 
of the order of 40%, which is caused by the suppression of the resolved contribution 
only. With the additional suppression of the DIRis term m me direct NLO correction, 
the ^-dependence of da/dQ 2 is reduced to approximately less than 20%, if we com- 
pare the maximal and the minimal value of da/dQ 2 in the considered £ range. The 
remaining ^-dependence is caused by the NLO corrections to the suppressed resolved 
cross section and the evolution of the virtual photon PDFs. How the compensation of the 
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FIGURE 1. Left: Photon factorization scale dependence of resolved and direct contributions to da / dQ 2 
together with their weighted sums for (i) suppression of the resolved cross section and for (ii) additional 
suppression of DIRis, using SaSID virtual photon PDFs 1 18]. Right: Q 2 -dependence of the dijet cross 
section for My = (full) and My = \\f T (dashed) and comparison with preliminary HI data using 
SaSID virtual photon PDFs 111 . 



My-dependence between the suppressed resolved contribution and the suppressed direct 
NLO term works in detail is exhibited by the dotted and dashed-dotted curves in Fig. 
[U (left). The suppressed resolved term increases and the suppressed direct NLO term 
decreases by approximately the same amount with increasing b ) . In addition we show 
also do/dQ 2 in the DIS theory, i.e. without subtraction of any ln<2 2 terms (dashed line). 
Of course, this cross section must be independent of £ . This prediction agrees very well 
with the experimental point, whereas the result for the subtracted and suppressed theory 
(full curve) lies slightly below. We notice, that for My = p\ the additional suppression 
of DIRis has only a small effect. It increases do /dQ 2 by 5% only. 

In order to get an idea about the My scale dependence of do/dQ 2 for the other Q 2 
bins we have computed this cross section for two choices of My, namely My = p^/4 
and My = Ap* T corresponding to the lowest and highest t, in Fig. [T] (left). The result for 
the do/dQ 2 is shown on the right side of Fig.[T] We see that the My-dependence in the 
considered range decreases with increasing Q 2 . This is to be expected since the resolved 
contribution diminishes with increasing Q 2 , so that the NLO corrections to the resolved 
cross section and the effect of the evolution of the photon PDF diminish as well. 



CONCLUSION 

In Summary, we described in this Report a new factorization scheme for diffractive 
production of jets in low-<2 2 deep inelastic scattering. By suppressing not only the 
resolved photon contribution, but also the unresummed logarithm as well as scheme- 



dependent finite terms in the NLO direct initial state correction, factorization scheme 
and scale invariance is restored up to higher order effects, while at the same time the 
cut-off invariance required in phase space slicing methods is preserved. 

For pedagogical reasons, we have chosen in this Report the kinematic region of finite, 
but low photon virtuality Q 2 , which exposes and regularizes a logarithmic virtual photon 
initial state singularity. We do, however, not rely on the finiteness of Q 2 , but rather 
separate suppressed and unsuppressed terms using the hard transverse momentum scale 
Pj, so that our scheme is equally valid for real photoproduction. 

The scheme- and scale invariance has been demonstrated numerically using the kine- 
matics of a recent HI analysis, differential in Q 2 and parton momentum fraction in the 
pomeron zp (not shown). Very good stability with respect to scheme- and scale varia- 
tions and good agreement with the experimental data has been found. 
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